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Abstract 

We consider random conformally invariant paths in the complex 
plane (SLEs). Using the Coulomb gas method in conformal field the- 
ory, we rederive the mixed multifractal exponents associated with both 
the harmonic measure and winding (rotation or monodromy) near such 
critical curves, previously obtained by quantum gravity methods. The 
results also extend to the general cases of harmonic measure moments 
and winding of multiple paths in a star configuration. 
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1 Introduction 



1.1 Historical perspective 

The subject of conformally invariant (CI) random curves in two dimensions 
has seen spectacular progress in recent years thanks to the invention of the 
Stochastic Loewner Evolution (SLE) (US [3]. This represents the crowning 
achievement of studies of 2D critical systems undertaken more than thirty- 
years ago. The first breakthrough came with the introduction of the so-called 
Coulomb-gas (CG) formalism. The critical properties of fundamental two- 
dimensional statistical models, like the O(N) and Potts models, could then 
find an analytic description within that formalism, which led to a profusion 
of exact results [U Ej . 

This was soon followed by the conformal invariance breakthrough that 
occurred in 1984 with the celebrated BPZ article [6]. It was followed by 
innumerable studies in conformal field theory (CFT), which became an es- 
sential source of applications to 2D statistical mechanics [7J [8]. That finally 
caught the attention of mathematicians [9], first through the peculiar cases 
of Cardy's formula for crossing probabilities in percolation [TU] and of inter- 
section properties of planar Brownian paths [TTj, resulting several years later 
in the advent of the SLE era [T2J ES [H] • 

Meanwhile, the deep relationship between the CG and CFT approaches 
was brought to light when the Coulomb gas representation appeared as an 
explicit model for a continuum of two-dimensional CFTs [TjJ [TB] . A Gaus- 
sian free field theory in a given domain, modified by a background charge 
2a (the "charge at infinity" ) that couples the field to the domain or bound- 
ary curvature, provides a concrete representation of abstract conformal field 
theories with central charge c = 1 — 24«q < 1 [T6]. Many of the critical 
properties of statistical models could then be obtained from the fusion of 
these approaches. 

Another breakthrough came in 1988 with the intrusion in statistical me- 
chanics of two-dimensional quantum gravity (QG) [17]. The famous KPZ 
relation between conformal weights in presence of a fluctuating metric and 
those in the Euclidean complex plane [TTJ [T8] could then be checked by ex- 
plicit calculations [T9J [2D1 121] • Though it is not the subject of this article, 
one cannot avoid mentioning in passing the enormous body of knowledge 
accumulated since then in the related field of random matrix theory. 

A decade later, just before the advent of SLE, it became clear, first from 
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the reinterpretation of an independent rigourous study in probability theory 
on intersection properties of planar Brownian paths [22] , that an underlying 
quantum gravity structure played an unifying role in conformal random ge- 
ometry in the complex plane [23]. In particular, revisiting QG allowed the 
prediction of fine geometrical, i.e, multifractal, properties of random critical 
curves. 

Those concerned the multifractal spectra associated with the moments of 
the harmonic measure, i.e., the electrostatic potential, near CI curves, which 
could be derived exactly within that approach for any value of the central 
charge c [21]. A generalization concerned the peculiar mixed multifractal 
spectrum that describes both the harmonic measure moments and the indef- 
inite winding or rotation (i.e., as logarithmic spirals) [25] of a CI curve about 
any of its points or of the Green lines of the potential [26J. Higher multi- 
fractality spectra were also introduced and calculated, concerning multiple 
moments of the harmonic measure and winding in various sectors of multiple 
random (SLE) paths in a star configuration [271 [28] . 

All these studies resorted to quantum gravity, by using a probabilistic rep- 
resentation of harmonic measure moments in terms of collections of Brownian 
paths, and by performing a "transmutation" of the latter paths into multiple 
(mutually-avoiding) SLEs, the rules of which are established within the QG 
formalism. The mixed harmonic-rotational spectrum was then obtained by 
blending this method with an earlier Coulomb gas study of the winding angle 
distributions of critical curves [2S] • 

In mathematics, the same multifractal harmonic measure or mixed spec- 
tra are the subject of present studies and can be obtained rigorously via the 
probabilistic SLE approach [30| I3T]. 

Recently, the same problem of the geometrical properties of critical curves 
was addressed in the physics literature via the Coulomb gas approach alone 
[32l [33] . That work in particular used the statistical equivalence between 
correlation functions of conformal operators in the complex plane C and 
correlation functions of a subset of these operators in presence of SLEs, i.e., 
in C cut by the latter [3H [33 ES]- Riemann uniformizing conformal maps 
were then devised to unfold the random paths onto the outside of some 
smooth domain, e.g., the unit disk. Standard CFT transformation rules of 
primary conformal operators were then the tool of choice. They allowed in 
particular to recover the multifractal exponents associated with harmonic 
measure moments near critical curves, originally obtained from quantum 
gravity. 
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This ingenious approach thus raised the interesting open question of how 
to generalize it to mixed harmonic and winding properties of random confor- 
mal curves, and recover the associated mixed multifractal spectra. The aim 
of the present article is to present such a generalization, still within the sole 
Coulomb gas formalism. 

Not surprisingly, the main ingredient is the use of arbitrary chiral primary 
operators, instead of the more familiar spinless ones. Their associated CFT 
transformation rules under arbitrary conformal maps will then yield informa- 
tion about the wild rotations that occur along Green lines near conformally 
invariant random paths. 

However, a technical subtlety will arise here, since one can no longer use 
only correlation functions of operators to extract geometrical information 
about windings of random paths. To study such windings, one has instead to 
resort to products of vertex operators, without statistical averaging, and study 
their monodromy properties at short distance. These encode the asymptotic 
geometrical rotation properties of the random paths and can be analyzed via 
operator product expansions, as will be shown below. 

1.2 Harmonic measure &; winding moments 
Definition 

Specifically, we consider a random conformal path (SLE) S and the harmonic 
measure u>(0, \z\) of a ball of radius \z\ centered at point G S, together 
with the (possible indefinite) rotation angle d(z) of the Green lines (i.e., 
electrostatic field lines) or, equivalently, of the equipotentials, when a point 
z tends to while avoiding S. We shall evaluate the asymptotic behavior for 
z — > of the mixed moment, 

u h {\z\) exp(-p&{z)) ~ \z\ £(h ' p) , (1) 

when one averages over configurations of the random path. The critical 
exponent x(h,p) depends on the two parameters h and p explicitly [26] . 

It is also interesting to consider in general multiple SLE paths with a star 
topology, and to define in a similar way the multiple mixed moments [25] : 

Ju^dzl) x exp(-p0(z)) ~ (2) 

where now co>j(|z|) is the harmonic measure in each sector i of the star, within 
a ball of radius \z\ centered at apex 0, while fl(z) is the common monodromy 
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angle describing the rotation for z — > of the Green lines between the arms 
of the random conformal star. Here again, the critical exponent x({hi},p) 
depends explicitly on all arbitrary exponents hi and rotation parameter p 
[28] , as we shall see in section [3j 

Riemann map 

A known way of computing these moments is to consider the conformal map 
w(z) which transforms the exterior of the random path S into the exterior of 
the unit disk D. The derivative of this map, w'(z), encodes all the relevant 
geometrical information. For 2; — s- 0, one has the equivalence ui(z) ~ |2«/(2)|, 
while the winding angle is given asymptotically by $(z) ~ — argw'(z). Thus 
the mixed moments above can be studied as well via the moments associated 
with the derivative of the conformal map 

\w\z)\ h exp[pargu/(2)] ~ \z\ x{h ' p \ (3) 

with now an obvious shift of the critical exponent x(h,p) := x(h,p) — h. 

In the case of a star, one takes a set of points 2, each in a separate sector, 
and with distances \zi\ ~ \z\ all scaling in the same way. Then the windings 
in each sector are equivalent so that for all i, arg 11/(2$) ~ — #(2), whence 

H\w'(z t )\ h >ex V {p ia rgw'(z t )} ~ \ z \<^\ (4) 

i 

where x({hi};p) := x({hi};p) - Y,i h i and P = J2iPi- 

To keep the formalism and technical notations to a minimum, and to avoid 
confusing the reader, we have simply adopted the notations of previous work 
by the Chicago group [35] . of which the present study can be considered as 
an extension. The reader is thus referred to their article which contains many 
relevant introductory details. The connection with our previous results and 
notations more familiar in quantum gravity will be recovered at the end of 
this article. 
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2 Derivation 



2.1 Star && operator products 

Vertex functions 



In complex Gaussian free field theory [371 EEL the basic objects considered 
here are "operator products", i.e., products of so-called vertex functions: 



P :-- 



\O a >. A >.(zi, Zi) 



i=l 



(5) 



where each vertex function or "operator" O a>& (z,z) is formally made of two 
respectively holomorphic and anti-holomorphic components 

O a , a {z,z) = V a {z)xV & {z) (6) 
V a (z) "= e^"^ 2 ) V & (z) ■= e'v^a^g) 

with Gaussian free-field correlators: 

{<j>{z)<j>{z')) = - ^g(z - z'), ($(zW)) = - log(« - z'), (<f>(zW)) = 0. (7) 

The holomorphic and the anti-holomorphic weights of the vertex operators 
are found in a standard way by applying the stress-energy tensor to them: 

h = h a = a(a — 2a ), h = h & = d(d — 2a ), (8) 

where 2a is the background charge acting in the Coulomb gas representation 
of the Gaussian free field theory. A given weight h thus corresponds to two 
possible charges: 



a h = a ± \]al + h, (9) 

an equation which applies separately to holomorphic and anti-holomorphic 
components. A vertex operator is spinless (meaning that h — h) if either 
a = a or a = 2a^ — a. In the sequel, since we are interested in the rota- 
tion of conformally invariant curves, the main tool is the consideration of 
chiral operators, as opposed to the usual case of spinless ones. Therefore we 
explicitly consider in (JSJ) operators O a ' i& '. with different conformal weights: 
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Star operator product 



To study the harmonic or rotation properties of a star configuration of n 
critical curves, i.e., SLEs, the main mathematical object is the star operator 
product [33] 



^n{{Zi,Zi}) :-- 



*0,n/3(°)li O °J,ajl 



2 = 1 



(10) 



J c 



where the operator ^0^/2(0) in ( TlQl) is, in the conformal field theory parlance, 
the operator corresponding to the existence of n critical curves originating at 
point z = in the plane. It represents the seed of a star configuration of n 
such curves, i.e., SLEs. The test operators O^a 1 - are a priori chiral, so that 
a- 7^ a- (or 7^ 2a — a-). 

We shall need the standard CG notations for holomorphic charges [38J: 



s )a- 



(11) 



where the basic charges are given in terms of the SLE parameter k g [0,8] 
by [33] 



H 



<K 2 

2a = + = — -=. 

1 \ K 



The charges a± satisfy the simple relations 



(12) 
(13) 



a± = a ± \Jal + 1, 

thus both correspond to a conformal weight h a± = 1, leading to the possibil- 
ity of using them as fundamental "screening charges" to build up the algebra 
of screening operators \W\ . 

The central charge of the CFT associated with the Gaussian free field 
theory, modified a la Feigin-Fuchs & Dotsenko-Fateev by the background 
charge 2a , is given by 



24ai 



3 (k 
2 



Ay 



K 



(14) 
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while the conformal weights corresponding to charges ct r ^ s are 

. „ . (rn - 4s) 2 - (ac - 4) 2 . . 

h r , s ■= a v (a r , s - 2a ) = — . (15) 

In this representation, the holomorphic (or anti-holomorphic) charge of 
the bulk curve-creating operator \&o,n/2 can be taken as ao,n/2 = «o — \ 0L - 
(with conjugate ao-n/2), and corresponds in the Coulomb gas formalism to 
a combination of electric and magnetic charges [531 EH]- The corresponding 
operator is spinless with conformal weight [39] 

W (16) 

Notice that operator products such as ([5]) or (fTUl) are understood as ob- 
jects to be inserted in correlation functions, in the same way as the formal 
holomorphic vs anti-holomorphic factorization (jfij) is meaningful only in such 
correlations (jSJ). Now, in the Coulomb gas formulation of conformal field 
theory, correlation functions are to be evaluated for a set of charges which 
overall respect electroneutrality, so that the sum of the latter always equals 
— 2ao, to compensate for the existence of the background charge 2a . So ul- 
timately, the set of charges introduced by the operator product P n near the 
origin, namely the set {a , n /2, (ca'i, = 1, ■ • - n} has to be compensated by 
other charges located at the observation points of the other remaining vertex 
functions. 

When dealing with "conformal blocks", namely treating separately the 
holomorphic and anti-holomorphic parts in correlation functions, electroneu- 
trality should apply separately to both sectors. Furthermore, if holomorphic 
(or anti-) electroneutrality does not apply in the arguments of the correlation 
functions of the original vertex operators, it is known that supplementary in- 
tegral screening operators, with screening charges a + and a_ ( TT2]) . can be 
introduced to extend the domain of definition of the theory. 

This caveat here is to explain why, when dealing with local operator 
products, such as that (jSJ) or (fT0|) . one does not have to enforce Coulomb 
gas neutrality rules locally, but only globally. We are here actually con- 
cerned only with the short distance and monodromy properties of the star 
operator product (I10p near the origin, which result from the analysis of the 
short-distance expansion of chiral operators O a >^' in presence of the spinless 
curve-creating one ^0,71/2; an d the CG formalism will be kept at a minimum. 
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It would be interesting to make a more complete study of the algebra of holo- 
morphic operators, including screening ones, for the present problem and in 
the context of SLE theory [161 HQ]- 



2.2 Star operator product expansion 
OPE for vertex functions 

Let us first consider an arbitrary product of vertex functions in the complex 
plane: 

p = n [°«,«i(*.*)]c = n [ v 4 z ^ x ^(*)] c ' 

where we have made explicit the formal factorization into holomorphic and 
anti-holomorphic vertex operators. 

When a subset V of I of points % G V is contracted towards a common 
point, like the origin here, i.e, Vz G V, z» — > 0, it is well-known that the 
limiting object is given by the so-called "operator product expansion" (OPE), 
with coefficients that are singular functions of the set of points V. 

In a Gaussian free-field theory, and for a product of vertex functions, i.e., 
complex exponentials of the field, these coefficients are precisely given by the 
average correlators of the set of vertex functions located at the contracted set 
V. The simplest and most convenient way to express this result is probably 
to use the normal ordering of operators : (• • • ) : such that [%Tl W2\ 

n p^mMc = <n ■ n \°<^^ c ■ ■ ^ 

i<=V ieV i£P 

In this notation, the vertex operators appearing inside : (• • • ) :, when inserted 
in any correlation function, are to be Wick-contracted only with operators 
located outside of contracting set V. 

In the complex plane, because of the form ([7]) of Gaussian averages 
(or 'Wick contractions", or also "free- field propagators"), the holomorphic 
and anti-holomorphic sector contributions to a correlation function decouple, 
hence factorize: 
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= (n[^)%(^)]>c 

= II {zi-z^izi-z^- (19) 

i,j€V,i<j 



This in turn gives the explicit form of OPE (1181) . 



Star OPE 

Applying the above OPE result to the star-operator product 



¥ n ({ Zi ,Zi}) 



i=l 



(20) 



when the set of points Zi, i — 1, ■ • ■ , n all contract to the origin z = 0, yields 

i<j 



Pn({^^})=^ ^ ^ o '" /2a ^ ao ^^ , 



z^z^^-z^ 



(21) 



x : F n ({zi,Zi}) : 



a result that will be our main tool in the analysis of the multifractal properties 
of a set of critical paths in a star configuration. 

2.3 Monodromy &; harmonic measure of a single path 
Original operator product 




Figure 1: Curve-creating operator \I>o,i(0) and chiral test operator O a ^(z, z), 
with monodromy of z around base point 0. 

Let us consider the standard case, where the star reduces to a single 
random curve passing through the origin 2 = 0, with n = 2 arms, i.e., a set 
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of two semi-infinite curves arising at the origin. The associated operator is 
^0,1 with holomorphic charge a ,i = a + /2. 

Let us place first one arbitrary "test" operator O a i^'{z, z) near the origin, 
so as to "measure" the harmonic measure moments and rotation on one side 
of the critical curve (FigflJ). The case of two test operators on both sides 
of the path will be treated later as the peculiar n = 2 case of the general 
n-star geometry. In the present geometrical situation, the relevant operator 
product is 

F(z,z) := [^ ,i(0)£W(^)] c , (22) 

and we are interested in its rotation, i.e., monodromy properties at short- 
distance (FigJT]). Its short- distance behavior when z — > is a simple case of 

(EU) 

F(Z, Z) = /«o,i°^2a ,ia' . ^ . ? (23) 

such that the SDE coefficient (123]) is explicitly 

¥(z, z) ~ z a+a z a+a = {zz) a+{a )/2 (y) . (24) 

Operator in path geometry 




Figure 2: Illustration of the statistical identity (125]) of the product of op- 
erators \I/o,i and O a) a with the latter one put in presence of the stochastic 
path. 

One first writes the identity in law of the operator product P ( 122]) in C, of 
the path creating vertex operator ^ ,i by the test operator, to the same test 
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operator but now in presence of the random path S originally represented by 
\l/o,i) hence to the same test operator in the complex plane C\«S slit by the 



curve: 



F(z,z) := [tt 0) i(0)<Va'OM)] 



(in law) 



(25) 



c 



c\s ■ 



This identity in law means that averaging within correlation functions the 
left-hand side over the complex Gaussian free field (GFF), or the right-hand 
side over the GFF in presence of S and over the configurations of the random 
path S yield the same result [3U [33]. [As discussed in [33], the precise 
boundary conditions on S for the geometrical random fields associated with 
the complex GFF may depend on the phase of the critical system ("dilute" 
for simple SLE paths, k < 4, or "dense" for k > 4), and are not specified 
here. In the holomorphic formalism, analytic continuation in k allows one to 
pass from one phase to the other.] 



Figure 3: Conformal map of the complement C\«S of the random path S in 
C to the exterior of the unit disk D. A mirror image of w(z) by inversion 
with respect to the unit circle <9D appears underneath the boundary. 

Conformal map 

The complex plane slit by S has the topology of the disk. One introduces 
the Riemann uniformizing map z — > w(z) that opens the slit S into the unit 
disk D centered at —i, so that w(0) = (FigJSJ) [321 E3]- This map naturally 
depends on the random path S. Under this conformal map, vertex operators 
transform like primary fields, whence [38J: 






[O a , t a>{z,z)] 



c\s 




C\B 



(26) 
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where the holomorphic and anti-holomorphic weights of operator 
respectively 



such that 



a (a 



a 



2a ), h' 



ha! 



a '(a' — 2a ), 



yi are 



(27) 



OLh' 



a Q ± 



h! 



Associating together the holomorphic and anti-holomorphic weights h! and 
h', we can write identically 



[O a ' A >(z, z)] 



c\s 



\w'(z)\ 



h'+h' 



dw \ («-*')/* 



0l i Al {w{z),w(z)) 



(28) 



where 



dw 



w'(z) 
w'(z) 



The final disk configuration in 



is a boundary con- 



O a > A i(w{z),w(z)) 

figuration [10], where an image charge appears at the conformal image of 
point w(z) with respect to the disk boundary, namely its image by inversion 
with respect to the unit circle. In the limit z —> 0, since w(0) = 0, the lat- 
ter image coincides with the image with respect to the tangent line, i.e., the 
complex conjugate w(z), so the same notation is kept here for simplicity. The 
vertex operator with anti-holomorphic charge a' in [0 '. n ']^\n then becomes 



an holomorphic vertex function of charge a' taken at image point w(z), so 
that the holomorphic and anti-holomorphic sectors get coupled |38j . 



SDE in the w plane 

When z — > in the original domain, w(z) — >• as well as w(z) — > 0, so the two 
image points pinch the unit circle at the origin in the w plane. The coefficient 
in the short-distance expansion (SDE) of the right-hand side operator in Eq. 
( 128]) is given by the Gaussian averaged correlation: 

(O a , A ,{w(z),^j)) = (w(z)-^{z)) 2 ^'. 
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Then the short-distance expansion of eq. (125]) is 

[O a , A ,{z,z)] cxs ~ \w'(z)\ h+h ^—(z)j (w(z) - w(z))™« . (29) 

After uniformization to the w plane, since the two image points w(z) and 
w(z) pinch the unit circle while staying either in the exterior or the interior of 
the disk, they cannot wind about the origin w(0) — indefinitely for \z\ — > 0, 
so arg u>(2:) remains bounded, e.g., 

argw(z) G [— 7r,+7r]. (30) 

We conclude that w(z) — w(z) ~ \z\ \w'(z) |, up to a non-winding (or non- 
monodromic) complex phase factor. This finally gives the SDE ([29]) 

/dm V")/ 2 
[O alA ,(z,z)] cxs ~ k| 2 ^V'(.)| h l^—(z)j , (31) 

where the overall (harmonic measure) derivative exponent h is defined as 

h:=ti + h' + 2a V, (32) 

Owing to fl27|) . it is also 

fa = (a 1 + a) (a + a — 2ct )- (33) 

It is the weight of an operator whose conformal charge ah = a® ± ^Jal + h 
is just a' + a' (or its conjugate): 



a' + a' = ah = «o ± y «o + ^- (34) 

Notice also that 

h' — h! — (a — a) (a + a — 2a ) . (35) 
Winding &c uniformizing map 

In the z plane, the curve S can wind or rotate indefinitely about any of its 
points, e.g. about the origin when z — > 0, so arg z — > ±oo. We are thus 
especially interested in the monodromy properties of the operator O a , „,(z, z). 
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By analyticity of the conformal map w onto the unit disk, w(z) ~ w'(z)z, 
so that a.igw(z) ~ arg w\z) + sxgz. We just have seen in (1501) that argio(z) 
remains bounded, so that asymptotically argz ~ — arg 11/(21) under the con- 
formal map. 

The winding angle ${z) of the Green lines of the random curve, asymp- 
totically close to z = on the curve, is given by 

•&(z) = argz = — axgw'(z) + 0(1). (36) 

Therefore one can also write asymptotically 

1 1 w'(z) 

$(z) = -Slogw'O) = - — [logiy'(z) -logw'(z)} = -— log==, (37) 

2i 2i w '(z) 

so that the exponential winding is 

where x here means equality within a non-winding (non-monodromic) phase 
factor. We therefore arrive at the expression for SDE ( 13T1) 



P*'AzMc\s ~ \z\ 2a ' a '\w'(z)\ h e-^ h '- h '^ . (39) 
Mixed moments 

Let us now return to the original OPE (|22j) . (123]) . We can rewrite (|23)» . (1231) 
as a complex scaling 

P(*,z) = [*0,l(0)O Q ',S^,2)] C ~ | z |OH-(o'+a') e <a + (a'-a')arg*_ ( 4Q ) 

Because of the identity in law (1251) . identifying (l4"0j) and the short-distance 
expansion (1391) in the transformed slit domain yields the equivalence 

This scaling equivalence is an SDE result, which can be interpreted as de- 
scribing typical statistical behavior, also expected to hold true in a weaker 
form after averaging over the configurations of the stochastic path S: 

-<\w'{z)\ h e- is ^ z V ~ \z\ x , (42) 
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where 

s := h'-h' + a+(a'-a') (43) 
x := a + (a' + d) - 2dd. (44) 

Let us now express exponent x := x(h,is) solely in terms of weight h and 
winding conjugate parameter s. Using (1331) (1541) and (1351) we have 

s = (d — a')(a+ + a' + a' — 2a ) = ( a ' _ a')(aq- + ah — 2a ) (45) 
x = a + a>h — lolol . (46) 

We can further write trivially 

a /-/ / / , -l\2 I I -i\2 2 I I -/\2 

4a a = [a + a) — (a — a) = a h — {a — a ) , 
to eliminate d — d between (l4"5l) and (1461) 



a' -a' = (47) 

a+ + ah — 2a 

1 1 s 2 

x{h,is) = a+a h --a 2 h + - — (48) 

z z (ct + + — za J 

To get the mixed harmonic measure-rotation exponents one has first 
to analytically continue s: s = —ip in ( 1421) so that 

^\w'(z)\ h e- p ^ z) y ~ |2r (h ' p) , (49) 

with now 

11 

a?(/i,p) = a+a h - -o^ - ■ — — -. (50) 

z 2 (a + + ah — 2a ) 

Finally, one has to choose the value of ah that vanishes with h, namely the 
"dilute phase" or simple SLE path one, where k < 4 and «o < (see (TT3"]) ) 



a ft = a + y«o + h. (51) 

The exponent x(h,p) above is identical to the scaling exponent obtained 
by quantum gravity in [26] for the hth power u h (z) e~ p ^^ of the harmonic 
measure uj with rotation conjugate parameter p (up to a natural shift by h, 
due to the local scaling uj(z) ~ for the harmonic measure w(z) in a 

ball of radius cf. Eqs. (JT| (j3J)). [See sectional] 
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2.4 Monodromy &; harmonic measure of multiple paths 

Operator product in star geometry 

We first use the identity in law of the star-operator product (120]) with the 
operator product of the test vertex operators in presence of the stochastic 
star S n in the complex plane: 

(52) 

C\5 n 



Pn(K*}) 



(in law) 



/ I ~? i Zi 




Figure 4: Conformal map of the complement C\c>3 of the random 3-star 
1S3 in C to the exterior of the unit disk D. Three mirror images of points 
w(zi),i = 1, 2, 3, appear by inversion with respect to the unit circle dD. 



Conformal map 



The plane slit by the star having the topology of the disk, a conformal map 
w(z) transforms the open set C\S n into the exterior of the unit disk D [33J. 
The vertex operator product P n fl52l) in presence of S n is transformed, ac- 
cording to primary operator rules, as: 



no 



c\s n 



H( W '( Zi )f (w'( Zi ) 



\0 <A (w(zi),w{zi)) 



(53) 
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Associating for each i the holomorphic and anti-holomorphic weights h! i and 
Kb we can write it identically 



dw 



(w(2i),«;(zi)) 



(54) 



c\i 



Short-distance expansion 

Now let all points converge to the origin in the original z plane, z$ — > 0, Vi = 
1, • • • , n. As in section (12.31) above, in the disk geometry, the operators on 
the right hand side of fj54jl are boundary operators, each of them made of a 
pair of holomorphic vertex operators taken at w(zi) and at its inverted mirror 
image with respect with the unit circle (fig. HJ). When Zi — > 0, each image 
coincides with the local complex conjugate w(zi) with respect to the local 
tangent to the unit circle, so the same notation is kept by a small abuse of 
notation. 

As usual, the original anti-holomorphic operators with charges ol\ now 
become holomorphic vertex operators located at points w(zi), and get cou- 
pled to the original holomorphic vertex operators with charges a[ located 
at points w(zi) |3Hj. Under the contraction Zj — > 0, Vz = 1, • • • , n, the SDE 
of the transformed operator product on the right hand side of (j54p there- 
fore scales as the Gaussian free-field average taken at all pairs of points and 
images [33] : 

(]JO <A (w(z i ),wJz7)))c\B = l[(w(z t )-wjz~)) 2 <^ (55) 

i i 

Among all these factors, the diagonal % terms give the short- distance be- 
havior. The % 7^ j terms are finite in the disk geometry, or give subleading 
contributions for configurations of the star in which a pair of points (zi,Zj) 
lie in the same sector, since the relative affixes w(zi) — w(zj), w(zi) — w(zj), 
w(zi) —w(zj) then tend to zero. The short-distance behavior of the operator 
product ( 1541) is thus 

„ - f dm XW-^)/ 2 

P„~]I \w'( Zl )\ h W ^(Zi) (w( Zi ) - w{ Zi )f<^. 
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When Zi — > 0, the transformed complex coordinate w(zi) and its image w(zj) 
pinch the unit circle, so that w(zi) —> w(zi) while the winding angle arg w(zi) 
stays bounded] whence we can simply take: w(zi) — w(zi) ~ up to 

a non-monodromic complex phase factor. This gives the final short- distance 
behavior of (154(1 



Pn ~ II (^(*) ) • (57) 

Hence the following exponents appear for the harmonic measure derivative 
terms : 

hi := h'i + h'i + 2a-^ = (a- + «•)(«■ + a • - 2a ), (58) 

which are the conformal weights of operators whose conformal charges simply 
result from the addition of holomorphic and anti-holomorphic charges : 



oihi = «o + y "o + hi = a'i + a'i. 
Windings & uniformizing map 

Let us now return to the SDE (l2Tj) in the original z plane 

p»(K m - n z?^z?^ H( Zl - Zj y^ (zi - Zj r^ (59) 

i i<j 



and, as in Eq. fpMj) . let us separate the complex modulus and argument 
parts, by rewriting it as : 



p»({*,*» ~ n { - 



Z; 



"o,»/a(«i-«i) 



^— ^) . (60) 

We are interested in the monodromy properties of expression (160]) when 
all points in their own sectors converge to the star apex z = while avoiding 
touching the star, which can wind indefinitely about its origin. Let us then 
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introduce a common scaling \z\ for all distances \zj\ to the star apex, as well 
as a common asymptotic winding angle d{z) for all arguments argZj: 

Zi = \z\e^Ci, axg* = 0{z) + argO, (61) 

where all moduli \Q\ and arguments aigQ remain bounded when \z\ — > 0. 
Since trivially 

(^) 1/2 = e^(0/0) 1/2 , (62) 
the short- distance expansion (1601) obeys the simple identity 

¥ n ({Zi,Zi}) = |^| 2Q o,n/2E 1 K+"i)+E i < j (2a-Qj+2a' i ^) 

x e «W[2a , n/a Ei(a' i -a' i )+2i; 4<j (a' i aJ-a{aJ)] 

x P n ({0,0})- (63) 

Let us now consider the geometrical setting in the transformed u>-plane. Since 
all arguments arg w(zi) stay bounded there, we have for Zi — > 0, as in Eq. 
( |36|) . aigw'(zi) = —argZi + 0(1), so that, as in (138]) 

Zi \ ^ 2 / du> 



w -l* w J ■ (64) 

Identity (1621) then yields immediately the common asymptotic behavior for 
all i 

( dw 



in terms of the unique rotation angle d(z) (JoTil , and up to non-monodromic 
phase factors. Applying this to the SDE result ( 1571) gives the asymptotic 
formula 

Pn „ ^(ZEi^e-wWEiW-'iiJJJ lu/te)!*. (66) 



Multiple moments and rotation 

Because of the identity in law (1521 . we can now identify expressions ( |63l) 
and (l66l) and get the fundamental scaling formula for the multiple harmonic 
measure factors and (indefinitely) rotating phase factor 

e -^W[E ! (^-^)+2ao,n/ 2 E ! K-^)+2E ! < J K«;-^^)] x JJ \ w '( Zi )\ hi 

i 

^ | z |2a ,n/2E,K+«D+E l<J (2^«;+2^a;)-E l 2^^_ (g-jA 
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We obtained this scaling relation from a short- distance operator product 
expansion, which yields the leading (typical) scaling behaviour of the product 
of harmonic measure factors and rotation phase factor. It can naturally be 
expected to also give the weaker result on the statistical average -<•••>- of 
the same quantity over the configurations of the multiple random paths of 
st ar i 

"41 W{zi)\ hi e-^ s y ~ \ z \<{hiY,is) ^ ( 68 ) 

i 

with the parametric representation in terms of holomorphic and anti-holomorphic 
charges 

hi := h' i + h' i + 2a' i a' i (69) 
h[ = a-(a- — 2a ); K = c^(c^ — 2a ) (70) 
s := - K) + 2a ,n/2 X)(aJ -aj) + 2 £(aja} - aja}) (71) 

x := 2a ,„/ 2 ^K + aO + X;( 2 «^ + 2 ^)-E 2 ^- ( 72 ) 

It remains to find the explicit expression of the scaling exponent x := x{{h{\\ is) 
in terms of the set of weights hi and of Fourier variable s conjugate to rotation 
Let us introduce the notations 

h a ' = olio! — 2ao) ; h a i = a (a — 2ao) 
a' := ^a'i, a' := J^aJ 

i i 

a h . = a'i + a'i 

a {h} := a' + a = ^(a- + a-) = ^Ja/*. 

i i 

A little bit of algebra first shows that s (I7ip can be written in the compact 
form 

S = k T,i oc'i + 2«0,n/2 Y a 'i ~ h T,t & 'i ~ 2a 0,n/2 «i 

i i 

= h a < - ha + 2a , n /2(a / - «') 
= (a — a)(a' + a — 2a + 2a ,n/2) 
= (a' - a')( Q! {/i} - 2a o + 2a ,n/ 2 )- 



(73) 
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One can also check that the exponent x (\72\i can be written in the compact 
form 

x = 2a 0i „ /2 ^ a hi + I a hi J - ^ a^. - 2a;' a' 

i \ i / i 

= 2a^ n/2 a [h} + - ^ a K ~ 2a '®'- ( 74 ) 

i 

Note that after this compaction of formulae, the expressions ( 1731) and ( 1741) 
are similar to expressions (|45j) and (146|) . We therefore eliminate a' — a' in 
the same way as above and arrive at a formula similar to (147|) and ( [48]) 

a' -a' = (75) 

ot{h} - 2a + 2a , n /2 



is) = 2a 0in/2 a w + ^a 2 {h} - ^ a 



1 s 2 

+ -7 v?, (76) 

2 (a w - 2a + 2a 0jn/2 y 



where we recall that a^y = a^, = a + a/«o + hi, hence exponent x 
has now an explicit form in terms of the set of weights {hi} and of s. 

It remains to analytically continue s into s = —ip to get the expectation 
of the multiple harmonic measure moments and Laplace transform of the 
rotation: 

-<Y[\w'(zi)\ h ^- p ^ ~ \z\ x{{h& ' p \ (77) 

i 

%({hi};p) = 2tt ,n/2 + \ a \h} ~ a hi 

1 p 



2 

-,2 



2 (a {h} - 2a + 2a , n / 2 )^ 



(78) 



CG formula (TTHj) for p = coincides naturally with the result found in ref. 
[33] . 
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3 Comparison to quantum gravity results 

In previous work [28] we introduced for SLE K the KPZ relation 

h = U K (A) := ^A (kA + 4 - k) 

between conformal weights A in a fluctuating metric with a conformal factor 
given by a Gaussian free field (two-dimensional "quantum gravity" (QG)) 
and conformal weights h in the complex plane. Its inverse reads 

A = U'\h) = ivi6Kfc + (K-4) 2 + i (\ _ 1 
Ik I \ k 

Using (TTH]) and (l5"Tj) we identify 



(79) 



The CG results above can be written as 



x({hi};p) = x({hi};0) 



2L {h} 



'{h} '■= ( a {h} - 2a + 2a ,n/2) 

1 o 



x({hi};0) = 2a , n/2 a {h} + ^a 2 {h} -^a 2 hi . 

i 

Using ()79p we find their equivalent quantum gravity expressions [28] 



»({^};0)-^ 



L{h} 

({hi};0) + J2 h i 



{h} 



2Z4 

2W W . 



1 /2 



2 



J {h} 



1 - 



1 [2 4 
- -n+ 1 

2 V At K 



(80) 



(81) 
(82) 

(83) 

(84) 



[As seen in (T4|) , the term £\ /ij simply comes from the passage from derivative 
moments to harmonic measure ones.l 
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Result (I50p for the multifractal exponent x(h,p) of a single-sided SLE 
simply corresponds to n = 2 and (hi, h 2 ) = (h, 0) in QG formulae ([8"Tl) - (l8"4"]) . 

The interpretation of these formulae is quite clear: 

• ()8ip is the Coulomb gas formula obtained a long time ago for the exponent 
governing the winding angle distribution of a star made of a number L^y of 
random paths (SLEs) [23 ESI M\ ; 

• L{h} (182]) represents, in a star topology, the effective number of SLEs that 
are exactly equivalent through QG to a collection of hi, i = 1, • • • , n, Brow- 
nian paths (which represent the set of powers hi of the harmonic measure), 
in addition to the n original SLEs [2"Tl |2"B"] ; 

• The scaling dimension (18"3|) is the result obtained through QG construction 
rules for (twice) the conformal weight in C of a random star made precisely 
of this effective number L{h} of SLEs, while (I84p is (twice) the weight of the 
original ra-SLE star without the auxiliary harmonic Brownian paths. (1831) 
is therefore also exactly (twice) the conformal weight h a — a(a — 2a ) cor- 
responding to the holomorphic charge a = ao,L {h} /2 of the curve-creating 
operator ^o t L {h) /2, which gives the final and rather elegant CG formula for 
flU: 

x({h,i}; 0) + ^2hi = 2a 0iL{h} / 2 (a 0>L{h} /2 - 2a ) - 2a ,„/2(ao,n/2 ~ 2a o)- (85) 

i 

Conclusion 

Results fl8TT) - fl84l) can be obtained by quantum gravity construction rules 
almost immediately, and have a natural interpretation in that formalism. As 
we saw, they can also be recovered under the form (I75|) in the fully developed 
Coulomb gas and conformal field theory approach, but though the latter is 
quite interesting, it also appears to be significantly more combersome. 

Its most important aspect is probably the renewed suggestion that a sys- 
tematic yet rigorous representation of the stochastic properties of SLEs via 
a Coulomb gas driven by a Gaussian free field must exist [43J, mimicking 
the (heuristic) physics CG approach's long-established predictive power. Let 
us finally mention that the mixed multifractal moments, introduced here 
to study the mixed harmonic measure-rotation spectrum of critical random 
paths, are amenable to an approach using only the Stochastic Loewner Evo- 
lution, that should allow one to establish rigorously the present results |30j . 
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While writing up this manuscript, we learned of related work by A. Be- 
likov, I. A. Gruzberg and I. Rushkin. 
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